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ABSTRACT 

A complete flux-limited sample of 118 ultraluminous infrared galaxies (ULIGs) has 
been identified from the IRASY&mt Source Catalog (FSC). The selection criteria were 
a 60 micron flux density greater than 1 Jy in a region of the sky 6 > —40°, \b\ > 30°. 
All sources were subsequently reprocessed using coadded IRAS maps in order to 
obtain the best available flux estimates in all four IRAS wavelength bands. The 
maximum observed infrared luminosity is Lir = 10^^'^'^ Lq, and the maximum 
redshift is z = 0.268. The luminosity function for ULIGs over the decade luminosity 
range Lij. = 10^^ — 10^^ Lq can be approximated by a power law oc L~^'^^ 

[Mpc~^mag^^]. In the local Universe {z 0.1), the space density of ULIGs appears to 
be comparable to or slightly larger than that of optically selected QSOs at comparable 
bolometric luminosities. A maximum likelihood test suggests strong evolution for our 
sample; assuming density evolution proportional to (1 + z)" we find a = 7.6 it 3.2. 
Examination of the two-point correlation function shows a barely significant level of 
clustering, ^(r) = 1.6 it 1.2, on size scales r ~ 22 Mpc. 

Subject headings: infrared: galaxies - galaxies: redshift - galaxies: luminosity function 
- galaxies: evolution - galaxies: clustering 



1. Introduction 

One of the important results of the Infrared Astronomical Satellite ^ (IRAS) all-sky survey 
was the discovery of a significant population of ultraluminous infrared galaxies (ULIGs: e.g. Soifer 
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et al. 1986, 1987; Sanders et al. 1987, 1988a) whose infrared luminositie^ (Lir > 10^^ Lq) are 
equivalent to the bolometric luminosities of opticahy selected quasars (e.g. Schmidt & Green 
1983). The origin of ULIGs and the source of their infrared power continue to be subjects of 
intense research. 

The IRAS 1-Jy survey of ULIGs was begun in early 1991 as an outgrowth of the original 
IRAS Bright Galaxy Survey (BGS: Soifer et al. 1986, 1987, 1989). The BGS, which included 
all extragalactic infrared sources brighter than 5.24 Jy at 60 ^m at high Galactic Latitude 
{\b\ ^ 30°), had resulted in the discovery of 10 ULIGs. The original intent of the 1-Jy survey was 
to increase this number of identified ULIGs by at least an order of magnitude in order to have a 
substantially larger sample of the nearest and brightest objects that could then be targeted for 
further ground-based study. 

The 1-Jy survey began by using the second version of the the IRAS Point Source Catalog 
(1988: hereafter PSC88). However, the PSC88 was abandoned in late 1992 following the release 
of the IRAS Faint Source Catalog - Version 2 (Moshir et al. 1992: hereafter FSC92). The 
FSC92 provided more accurate fluxes^] and positions than the PSC88 and is now considered the 
standard reference catalog for IRAS sources. Section 2 provides a complete description of the 
selection criteria used in the 1-Jy survey and describes how candidates were selected from the 
FSC92. Section 3 describes the redshift survey that was subsequently carried out to identify 
ULIGs. The final 1-Jy sample (including finder charts for each source) is presented at the end 
of Section 3. Section 4 derives the luminosity function for ULIGs, while Section 5 discusses the 
general properties of ULIGs and compares the luminosity function of ULIGs with other classes of 
extragalactic objects. 

More detailed analyses of individual sources in the 1-Jy sample are being presented in 
companion papers. Optical spectra of ULIGs that were obtained by Kim (1995) are analyzed in 
Kim, Veilleux & Sanders (1997: Paper II). More recently obtained optical spectra for the entire 
1-Jy sample are analyzed in Veilleux, Kim, &: Sanders (1998). Near-infrared spectra for ULIGs 
that were obtained with the United Kingdom Infrared Telescope (UKIRT) on Mauna Kea are 
presented in Veilleux, Sanders, &: Kim (1997; 1998). Optical and near-infrared images that were 
obtained by Kim (1995) for the majority of the ULIGs are analyzed in Kim, Sanders, &: Mazzarella 
(1998: Paper III). 



''l/ir = L{8 — lOOO/xm). Throughout this paper we assume Ho ~ 75 km s ^Mpc ^ ,qo — 0.0 

^The FSC92 incorporates the final calibration scale (c.f. Pass-3 calibration) that removed inconsistencies in the 
previous point source fluxes computed prior to 1991 versus those computed using the more accurate IRAS data 
products such as ADDSCAN/SCANPI (Helou et al. 1988). Pass-3 calibration had the largest impact on the IRAS 
60 fim and 100 ^m bands where typical flux changes, compared to values in the PSC88, were ~ 5-10 % with a maximum 
change of ~ 15 % for the sources in our sample. 
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2. The IRAS FSC Survey 

2.1. Selection Criteria 

Four selection criteria were used to extract candidate ULIG sources from the FSC92. Criteria 
1-3 follow the general selection methods adopted earlier for the IRAS BGS. An important fourth 
criteria was added following the detailed analysis of the properties of IRAS galaxies in the BGS 
by Soifer & Neugebauer (1991). The four criteria are: 

(1) All sources were required to have /eo > 1 Jy^, with a "high" or "moderate" flux quality index 
(Moshir et al. 1992) This limit is ~4c7 above the typical threshold at 60 /im for inclusion in the 
FSC. 

(2) As in the original BGS, all sources were required to he at high Galactic latitude (\h\ > 30°), 
and to be easily accessible for redshift measurements (in this case, 5 > —40° to allow for observing 
from Mauna Kea). The Galactic latitude restriction was imposed to avoid serious contamination 
by infrared emission from the Galactic plane. The total sky coverage represented by the current 
survey, excluding small regions of the sky not surveyed by IRAS or with less than two IRAS visits 
(see IRAS Catalogs and Atlases Explanatory Supplement) is ~ 16,300 sq. deg. Compared to the 
original BGS, the current survey area is approximately 10% larger due primarily to the lower 
declination limit used at Mauna Kea. 

(3) As in the original BGS, all sources were required to have /eo > fi2 in order to exclude 
infrared- bright stars which are known to peak around 12 /im (e.g. Cohen et al. 1987). This criteria 
has been shown to be a very efficient means of selecting for galaxies while weeding out most 
Galactic objects. Our experience with the BGS showed that the remaining sources of confusion 
are from relatively nearby planetary nebulae, heavily reddened stars, and reflection nebulae. All 
galaxies in the BGS have /go > fi2, and this requirement is particularly robust in the case of 
ULIGs where /60//12 is typically > 20 (see Soifer & Neugebauer 1991). 

(4) All sources were also required to have "warm" /eo/Zioo colors: log (/eo/Zioo) > —0.3. Soifer 
&Z, Neugebauer (1991) showed that for galaxies in the BGS, the mean /eo/Zioo flux ratio increases 
systematically with increasing infrared luminosity, and at the highest observed luminosities (i.e. 
for the 10 ULIGs in the BGS) the mean ratio is —0.11 with a range of —0.2 to +0.13. Subsequent 
analysis of an additional 25 ULIGs from both the 2-Jy Survey of Strauss et al. (1992) and the 
BGS-Part II (Sanders et al. 1995) revealed nearly the same mean and range of far-infrared colors. 
Therefore, it was decided that the criterion of log (/eo/Zioo) > —0.3 would not exclude legitimate 
ULIGs, and that it could indeed be used to exclude a large fraction of the more numerous 
lower-luminosity galaxies, the vast majority of which have cooler /eo/Zioo colors. 



® The quantities /12, /25, /so, and fioo represent the IRAS&ux densities in Jy at 12 pm, 25 /im, 60 fim, and 100 nm 
respectively. 
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2.2. Selection of ULIG Candidates 

Our search of the FSC92 using the four criteria listed above produced a list of 1630 candidates 
from a total of ~ 175,000 catalog entries. Two steps were then used to further reduce this number; 
first by cross-correlating our list with catalogs of known Galactic and extragalactic sources, and 
second by using a size criteria to eliminate relatively nearby lower luminosity galaxies. 

All 1630 IRAS sources were cross-correlated with previously published astronomical catalogs 
by using the catalogs available in the NASA/IPAC Extragalactic Database (NED) and the 
Strassbourg Data Base (SIMBAD III). This resulted in the exclusion of 576 sources that were 
either Galactic objects (mostly highly reddened stars, planetary nebulae, and reflection nebulae), 
or nearby lower luminosity galaxies. The NED search also picked up 28 previously identified 
ULIGs - 10 from the final version of the BGS (Soifer et al. 1989), and an additional 18 ULIGs 
from various published lists of IRAS galaxies. Another 20 ULIGs were then identified from a 
preprint version of the /eo > 1.936 Jy survey eventually published in Strauss et al. (1992). Thus, 
the net result of the cross-correlations was a list of 48 previously identified ULIGs and a reduction 
in the number of remaining unidentified ULIG candidates from 1630 to 1006. 

Palomar Optical Sky Survey (POSS) overlays were produced for the remaining 1006 
unidentified candidates, and the optical images were examined using the POSS prints. The 
majority of these objects appeared to be spiral or irregular galaxies as judged by their appearance 
on the POSS E-prints. Most of these galaxies had optical diameters that were much larger 
than any of the nearest known ULIGs, and were therefore assumed to be relatively nearby, 
less-luminous, infrared galaxies. 

The following simple procedure was developed in order to exclude nearby, lower-luminosity 
galaxies. An optical diameter twice that of Arp 220 (at z = 0.0185, the nearest ULIG by a factor 
of ~ 2) as measured from the POSS E-plate, was used to set a minimum size above which any of 
the remaining 1006 sources was excluded from the final list of ULIG candidates. This procedure 
eliminated 837 sources [|. However, 18 sources were put back in the unidentified category due to 
the fact that there were two extended objects within the IRAS position uncertainty ellipse, and 
one of these objects was smaller than our template diameter. This resulted in a final list of 187 
unidentified candidate ULIGs. 

Overlays were also produced for the 48 previously identified ULIGs in order to verify that the 
more accurate positions from the FSC92 were still spatially consistent with the source, all of which 
had previously been identified using the PSC88. All sources except for one were found to still lie 
well within the ±3 a FSC92 position uncertainty ellipse. That one source, IRAS F13442-I-2321 
(Leech et al. 1989) was found to be centered on a low luminosity spiral galaxy rather than the 
previously identified higher redshift galaxy which was more than 5 a away from the IRAS FSC92 



^ All of these sources were later proven to be lower luminosity galaxies by cross-checking with the published 
redshift surveys of Strauss et al. (1992) and Fisher et al. (1995). 
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position. This source was dropped from our list. 

3. Identification of ULIGs 

The following procedures were used to determine the final list of ULIGs: (1) Optical spectra 
were obtained of all 187 previously unidentified candidates. In cases were more than one obvious 
source was found within the IRAS position uncertainty ellipse spectra were obtained of each 
object. (2) luminosities were computed using measured redshifts from step 1 and the fluxes in 
all four IRAS bands. (3) optical spectra at higher spectral resolution were obtained of all ULIGs 
identified from step 2 and all 47 ULIGs identified from the literature (as part of a follow-up 
program to obtain line-ratio diagnostic information for all ULIGs). The final list of 1-Jy ULIGs is 
given in Table 1. 

3.1. Redshift Measurements 

Redshifts were initially obtained for all 187 candidate ULIGs and the 47 previously identified 
ULIGs using the University of Hawaii 2.2-m telescope on Mauna Kca (Kim 1995). Optical spectra 
were obtained during five separate observing sessions between 1991, August and 1993, March using 
the Faint Object Spectrograph (FOS) at the f/10 Cassegrain focus of the UH 2.2-m telescope. 
Wavelength coverage was typically A = 5, 000 — 9, 000 A at a resolution of ~ 5 A. Exposure times 
varied from 180 sec to 1800 sec depending on the optical brightness of the source. A Zeiss 
135R camera and long slit (2.5" x 275") were used for all observations except during February 
1993 when a Nikon 180R camera was used. All observations were obtained with either a Tck 
1024x1024 or Tek 2048x2048 CCD. The seeing was typically ~0.8"- 1.4". Data reduction was 
performed using the TWODSPEC package in IRAF. All spectra were de-biased by fitting a high 
order polynomial to the bias column, and flat-fielded with a normalized dome-flat. Wavelength 
calibrations were performed using an Fe-Ar arc lamp, and flux calibration was done by observing 
spectroscopic standard stars. The redshifts were determined from Gaussian fits of Ha, [Nil] and 
[SII] emission lines supplemented on occasion with fits of MgIb and NalD absorption lines. The 
typical uncertainty in the measured velocities is estimated to be ~ ±70 km s^-*^ (l^"")- 

Out of the 187 previously unidentified candidate fields, 175 were found to contain emission 
line galaxies. Of the remaining 13 objects, 4 were stars, and 9 turned out to be empty fields which 
most likely are infrared cirrus or small Galactic dark clouds. 

Six of the 175 fields were found to contain two equally probable galaxies within the ±3 a 
position uncertainty ellipse. In each of the six cases both galaxies had nearly identical redshifts 
{Az < 0.002), and high resolution imaging (Paper III) shows that all six appear to be strong 
tidally interacting systems; therefore for these six fields it was decided to attribute the IRAS 
source to the pair and to continue treating the source as a single entry in our list of ULIGs. 



- 6 - 



3.2. Infrared Luminosities 

Using our redshifts for each of the 175 newly identified extragalactic sources and published 
redshifts for the 47 previously identified ULIGs, an infrared luminosity, Lir, was calculated using 
the expression 

Lir = ^TTDj X Fi, (1) 

where is the luminosity distance for a Qq = cosmology (Weinberg 1972) 

De = -^{z{l + 0.5z)) (2) 

and the infrared fiux, Fjr, in units of [10^^^ W m~^] is defined by 

Fir = 1.8 X {13.48 X /i2 + 5.16 x /ss + 2.58 x /go + 1.0 x /loo} Q. (3) 

The IRAS fiuxes were taken directly from the FSC92. A substantial fraction of the sources only 
had upper limits listed for /12, and a smaller fraction had upper limits listed for 725. 

In an attempt to provide the most accurate luminosities possible, we attempted to 
remove 12 //m and 25 /im upper limits by manually reprocessing the IRAS data using the 
ADDSCAN/SCANPI procedure (Helou et al. 1988) and visually inspecting the output. For 
~ 50 % of the objects with 25 /xm upper limits, and for ~ 20 % of the objects with 12 ^um upper 
limits, we were able to find a signal at the 2-3 a level; for these cases it was decided to use the 
new values rather than continue to use the upper limits in the FSC92. For those upper limits that 
could not be removed using ADDSCAN/SCANPI we decided to use a template SED (constructed 
from a subsample of nearby ULIGs) to extrapolate the 25 /xm and 12 fxni fiuxes. The template 
used was fu = 0.20 /25, and /25 = 0.14/60 for objects with Lir < 10^^'^ Lq and /25 = 0.23 /60 
for objects with Lir > lO^'^'^ Lq (see §5.1). Given these procedures for removing upper limits at 
12 //m and 25 /xm, and the fact that even with the upper limits the sum of the 12 /im and 25 fim 
contribution to Lir is typically < 10 %, we estimate the final uncertainty in Lir due to the 12 jim 
and 25 fim flux densities to be a few percent, i.e. similar to the measurement uncertainty in the 
60 /xm and 100 fxm flux densities. 

Seventy-one of the 175 galaxies proved to have infrared luminosities large enough to be 
classified as ULIGs. The remaining 98 objects were lower luminosity emission line galaxies. 



*The factor of 1.8 in equation (3) was found by Perault (1987) to provide a good approximation (±5%) to the 
ratio of the 8-1000 fim flux (= Fir) and the sum of the fiux in all four IRAS bands (quantity in brackets), for thermal 
dust emission in the temperature range ~ 25-300 K and for dust emissivity laws e oc X~" for a in the range -0.5 to 
-2.0. Since many of the high-luminosity galaxies emit a significant portion of their infrared luminosity shortward of 
40 /xm, this expression provides a significantly better determination of the total infrared flux than the more commonly 
used Ffiy determined by fitting a single temperature dust model to the 60 /xm and 100 /xm fluxes (cf. Appendix B of 
Cataloged Galaxies and Quasars Detected in the IRAS Survey 1985). 
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3.3. Final Adjustments 

As part of a parallel program to obtain broader wavelength coverage plus higher spectral 
resolution data for all of the ULIGs in the 1-Jy sample (Veilleux, Kim, & Sanders 1998), we were 
eventually able to obtain our own redshifts for all of the 47 previously identified ULIGs while also 
being able to confirm our redshifts for the 71 previously unidentified ULIGs. For all 71 newly 
identified ULIGs. and for most of the objects taken from the literature, our new redshifts were 
in excellent agreement (Az < .001) with previous data. Notable differences {Az > .002) were 
found only for six objects, all of which were from the subsample of 47 objects taken from the 
literature. However, none of the 47 previously identified ULIGs fell out of the ULIG sample when 
their luminosities were recomputed using our measured redshifts and the new FSC92 fluxes. A 
comparison of our final redshift measurements with rcdshift measurements determined by others 
is given in Table 2. The final 1-Jy sample contained 71 + 47 = 118 ULIGs. 

3.4. The Complete IRAS FSC 1-Jy Sample 

Table 1 lists source parameters for all 118 ULIGs in the final list of objects in the 1-Jy 
survey. Column entries are described in the footnotes. All source parameters in columns 1-12 
are from the FSC92, except for the estimated 12 //m and 25 //m fluxes (FQ = a in column 12). 
In addition to the measured redshift (column 13) and computed Lir (column 15) we have also 
provided the quantity Lfij. (column 14) for comparison with previously published IRAS galaxy 
samples which used this quantity as the primary measure of "far-infrared" luminosity. A more 
complete description of the various methods used in the past to represent the IRAS luminosities 
is given in Sanders &: Mirabel (1996). 

3.4- 1- Finder Charts 

To aid in the correct identification of the actual IRAS source for each of the objects in Table 
1, we have provided small finder charts in Figure 2 showing the IRAS position uncertainty ellipse 
properly registered on the red filter of the Palomar Digital Sky Survey with a marker to identify 
the ULIG. The distribution of offsets between the actual source and the IRAS position listed in 
the FSC92 (columns 2-5 of Table 1) is typical of the ~ 10-15" (lo-) uncertainty quoted for the 
FSC. All of the sources in the 1-Jy survey are located within the ±3a uncertainty ellipse. 

4. The Luminosity Function 

The rcdshift and luminosity distributions for the 118 ULIGs listed in Table 1 are plotted 
in Figure 3. The mean and median luminosity values for the sample are (i^ir) = lO^^'^^^'^'^^ Lq 
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and Lir = 10^^-^^ Lq respectively. The mean and median redshift values for the sample 
arc {z) = 0.144 it 0.043 and z = 0.145 respectively. The maximum observed luminosity is 
(-^ir)max — ]^q12.90 ]_,q^ and the maximum observed redshift is ^imax — 

0.268. 

The infrared luminosity function of ULIGs was calculated using the following expressions 
(Schmidt 1968) 



and 

'47r\ / 1 



where Q, is the solid angle of the survey region, and is the maximum volume within which 
the object would still be included in the survey. V and V^ax were calculated by determining the 
luminosity distance according to the following procedure. The solar motion with respect to the 
Local Group — 300 sin £ cos 6 [km s~^] (de Vaucouleurs, de Vaucouleurs, &; Corwin 1976), where £ 
and h represent Galactic longitude and latitude respectively — was subtracted and the resulting 
velocity was then converted to a luminosity distance using equation (2). A correction for the Virgo 
flow was not included since most of our objects are at distances far beyond the influence of the 
Virgo cluster. The monochromatic flux density can then be written as (Weinberg 1972) 

- 4vrD// (1 + z) 

where Ly(^ij^z) is the intrinsic luminosity at frequency v{\ + z) and is the luminosity distance. 
From equation (6), we can then calculate a maximum luminosity distance Z)^(max), which is 
related to the limiting flux density Fmin by 

D,{ms^) = DAJ—-^-^^) (7) 

\ /min -L I -Zmax / 

where the i^'-correction term ^ was calculated by a polynomial fit to the spectral energy 
distribution between 25 /im and 100 /xm. The typical value of ^ in our sample was found to be 
0.92-0.94. The luminosity function, <1> [Mpc^'^ M/2^^], for the 1-Jy sample is tabulated in Table 3. 



5. Discussion 

The luminosity function for ULIGs provides a basis for discussing the relative importance 
of these sources in relation to other populations of extragalactic objects. However, before 
discussing the luminosity function of ULIGs it is important to understand just how dominant 
is the far-infrared emission from these objects, and why they have until now mostly escaped 
identification in optical surveys. 
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5.1. Spectral Energy Distributions of ULIGs 

The IRAS fluxes from Table 1 have been combined with optical and near-infrared fluxes 
that have recently been obtained for a substantial fraction of the 1-Jy sample (Kim 1995, Kim et 
al. 1998), to compute the mean SEDs of ULIGs as shown in Figure 4. The mean SED for the 
entire sample, (as well as the SED for each ULIG in the sample), peaks in the far- infrared near 
60 //m. For all of the ULIGs in our sample, the total infrared luminosity, L^^, appears to be by 
far the dominant fraction of the bolometric luminosity. Eventhough the relative optical/infrared 
luminosity ratio increases for objects with Ljr > 10^^'^ Lq (see Figure 4), Lir/Lopt+ir is still >0.9 
for all but three of the ULIGs in the 1-Jy sample. Given that the amount of energy outside the 
infrared band is typically < 10 % of Lji- in all but a very few of the 1-Jy ULIGs, we have used Lir 
as a good approximation to the bolometric luminosity for the purpose of comparison with Lboi 
calculated for other extragalactic objects. 

Compared to modest infrared galaxies like our own (Lir ~ 10^° L0, Lir ~ Lopt), ULIGs have 
larger infrared luminosities by at least a factor of 100 while their optical luminosities have only 
increased by factors of ~ 2~5. The mean optical and near-infrared luminosities of ULIGs are in fact 
~ 2.5 L* in the K-band 0, ~ 2.7 L* in the r-band, and ~ 2 L* in the B-band (see Sanders & Mirabel 
1996). Given the fact that the nearest ULIGs (i.e. those in the RBGS) have mean redshifts ~0.04, 
none of these objects are included in optical catalogs of bright galaxies (e.g. NGC), and only four 
objects in the 1-Jy sample are found in the fainter UGC and Zwicky catalogs. Thus the IRAS 
ULIGs represent an independent sample when compared against the well-known optical catalogs 
of galaxies. 

Further discussion of the SEDs of ULIGs, including a comparison of SED properties with 
object morphology and spectral type, will be given in Papers II and III. Here it is simply worth 
noting that the only obvious correlation of the shape of the infrared SED with galaxy properties is 
an apparent increase in the /25//6O ratio with increasing infrared luminosity. A significant fraction 
(~ 1/4) of the ULIGs with Ljr > 10^^'^ Lq have /25//6O > 0.2, similar to values that are typically 
found in infrared selected AGN (e.g. Miley et al. 1984; de Grijp et al. 1985; Golombeck, Miley 
& Neugebauer 1988). Paper II discusses the possibility that the increase in the /25//6O ratio is 
in fact due to an increasing proportion of AGN (including a substantial fraction of Seyfert Is) at 
the highest luminosities in our ULIG sample, a conclusion that is supported by new near-infrared 
spectra of a subsample of 25 ULIGs chosen from the 1-Jy survey (Veilleux, Sanders, &: Kim 1997). 



^ For an L* galaxy Mg = —19.7 mag (Schechter 1976), M* — —20.5 mag assuming a typical B — r — 0.75, and 
Mk- = -24.2 mag (Mobasher, Sharpies, & Ellis 1993). 



- 10 - 



5.2. Comparison of the Luminosity Function of ULIGs with that of Other 

Extragalactic Objects 

Soifer et al. (1987) presented the first meaningful comparison of the luminosity function 
of infrared galaxies with other classes of optically selected extragalactic objects by calculating 
a bolometric luminosity function for both infrared and optical samples. Figure 5, which is an 
updated adaptation of a figure from Soifer et al. (1987), illustrates the differences and similarities 
between the more well-known optical classes of galaxies and the infrared galaxies revealed by IRAS, 
and shows how infrared galaxies become an increasingly dominant population with increasing 
bolometric luminosity. Except for the new 1-Jy sample, the comparisons illustrated in Figure 5 
only reveal what is found for the "local Universe" (i.e. z ^ 0.1) in that the redshift surveys of 
the optical samples as well as the RBGS typically do not extend much above a mean redshift of 
~ 0.05. The one exception is for QSOs, but this has been dealt with in Figure 5 by showing the 
local luminosity function for QSOs, ^{Mb,z = 0), as derived by Schmidt & Green (1983). 

The luminosity function for normal galaxies (i.e. the Schechter function) has its inflection 
point near Lir ~ 10^'''^ Lq (i.e. L*), and falls exponentially at Li^ > 10^^ Lq. The luminosity 
function for optically selected starbursts is remarkably close to the power-law defined by the 
infrared galaxies at Lij. < 10^^ Lq, but starts to fall precipitously at higher luminosities, with 
few identified objects at Lji- > 10^^'^ Lq. On the other hand, optically selected Seyfert galaxies 
(primarily type-2 Seyferts) are less numerous than both infrared galaxies and optically selected 
starbursts at Ljr < 10^^ Lq, but are approximately equal to the space density of infrared galaxies 
over the decade luminosity range Li, = W^^-^^ Lq. At Lboi > Lq, only QSOs (i.e. type-1 
Seyferts with Mr < -23, for Ho = 75kms"^ Mpc~^) remain in the optical samples, and figure 5 
shows that the space density of PG QSOs in the local Universe is ~ 1.5 times lower |^ than that 
of ULIGs in the RBGS. However, the relative space densities of optically selected QSOs may be 
equal to or even slightly exceed that of ULIGs if one accepts the most recent results of Wisotzki 
et al. (1998; see also Krishna & Biermann 1998), who claim to find larger numbers of QSOs in the 
local Universe (compared to the number density of PG QSOs) from their new flux-limited sample 
of bright low-redshift QSO and Seyfert 1 galaxies drawn from the Hamburg-ESO survey. 



5.3. Evolution 

Figure 5 also shows the space density of ULIGs from the 1-Jy sample. The space density of 
ULIGs in the 1-Jy sample appears to be ~ 2 times larger than that computed for ULIGs in the 



^"The original comparison of ULIGs and QSOs (Soifer et al. 1987) found that the space density of QSOs in the 
local Universe was ~ 3 times smaller than that of ULIGs in the BGS. This difference was reduced to a factor of ~ 2 
by Sanders et al. (1989) after better accounting of the infrared emission from QSOs, and finally to the factor of 1.5 
reported here due to using a better approximation to the far-UV and soft X-Ray component of QSOs as given in 
Elvis et al. (1994). 
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IRAS RBGS. The 1-Jy ULIG survey samples a larger volume of space (-S^max OS with, a mean 
ULIG redshift of z ~ 0.15) than the RBGS (zmax ~ 0.08 with a mean redshift of z ~ 0.04 for 
ULIGs in the RBGS). We show below that this increase in the space density of ULIGs is consistent 
with strong evolution in the ULIG population. 

To quantify the evidence for evolution in the luminosity function of ULIGs we used the 
maximum likelihood test of Saunders et al. (1990) including refinements by Fisher et al. (1992). 
We also assumed pure density evolution partly because our sample is not cosmologically distant 
enough to differentiate between different evolutionary models (e.g. pure density or pure luminosity 
evolution), but also because the most recent results from AGN surveys appear to favor pure 
density evolution over pure luminosity evolution models (e.g. Hassinger 1998; Kohler & Wisotzki 
1997). The probability pi that a galaxy i of given luminosity is located within a comoving volume 
dV centered on its observed position with an isotropic comoving density distribution N(r) is 



If we assume that the comoving density distribution N{r) has the form N{r) = A^(0)(1 + z)", 
where A^(0) is the local space density of ULIGs, then the overall probability of measuring the 
ensemble of a total of n galaxies each with position is 



For an assumed cosmological model with given object redshifts and luminosities, a maximum 
likelihood estimate for a can be found by maximizing the probability P with respect to a. In this 
calculation, the maximum distance rmax was calculated using equation (7) where a polynomial fit 
to the spectral energy distribution is used to determine the _R'-correction W'. For = 0.0 and Hq 
= 75 (km s~^ Mpc~^), the formal maximum likelihood estimate for our sample was a = 7.6 ± 3.2 
This suggests strong evolution in the ULIG population. However, it should be noted that the 
redshift range over which this estimate applies is relatively small {z ^ 0.27), and the corresponding 
uncertainty in the exponent is relatively large. 

Our value of q = 7.6 ± 3.2 assuming a pure density evolution model can be compared with 
previous estimates of evolution that have been determined from other redshift surveys of IRAS 
galaxies. (However, we note that no previous analyses have focused only on ULIGs, but instead 
have considered all infrared galaxies regardless of luminosity.) The most relevant measure of 
evolution for our purposes comes from a recent analysis of the 'FSS-z' sample (Oliver et al. 1996) 
of sources with /go > 0.2 Jy from the IRAS Faint Source Data Base. Oliver et al. (1995) used 
~ 1200 sources over the redshift range 0.02 < z < 0.3 for which redshifts were available, and found 
a = 5.6 lb 1.6. The redshift range is nearly identical to that of our 1-Jy sample of ULIGs, but the 
FSS-2; sample, of course, includes all sources regardless of luminosity. Oliver et al. (1995) also 
provide a rather thorough accounting of several other determinations of a in the literature, all of 
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which have been made using brighter 60 /xm samples, and all of which, when properly corrected for 

Malmquist bias, find lower values still for a, including a suggestion of no evolution, e.g. a = 2 it 3, 
for the /go > 1-2 Jy sample (Fisher et al. 1992) where no lower redshift cutoff was imposed. This 
latter case is clearly affected by the known local overdensity of infrared sources (e.g. Hacking, 
Condon, k. Houck 1987), which when included in the analysis would mimic negative evolution 
(i.e. the value of a increases to 4.2ib2.3 for this same sample just by imposing a lower redshift 
cutoff, z > 0.02). Although these various measures of evolution at first appear inconsistent, they 
are in fact consistent with the idea that strong evolution is only found for the most luminous 
IRAS galaxies, and then only for rcdshifts beyond the influence of 'local' large-scale structure (e.g. 
z ^ 0.05). Similar conclusions can also be drawn from previous analyses of IRAS extragalactic 
source counts (Hacking et al. 1987, Lonsdale & Hacking 1989, Lonsdale et al. 1990, Gregorich et 
al. 1995) that show evidence for strong evolution only at relatively low flux levels (/eo < 1 Jy)- 

It is intriguing to note that strong evolution has also been found for optically selected QSOs. 
Schmidt Sz. Green (1983) find clear evidence that the strength of evolution depends on luminosity, 
with little or no evidence for evolution in objects less luminous than QSOs (i.e. Mb > —23 

corresponding to Lboi < 10^^ Lq) with evolution rapidly becoming stronger for more luminous 
systems. Although Schmidt &: Green (1983) discussed their sample in terms of pure luminosity 
evolution, an examination of their data assuming pure density evolution shows that the luminosity 
function evolves as (1 + z)~^~^ for QSOs with redshifts and luminosities similar to what is 
sampled by the 1-Jy survey of ULIGs (i.e. for QSOs with Mb = —23 to —25 and redshifts in 
the range z = — 0.3). Unfortunately, there are insufficient data on higher luminosity ULIGs 
(i.e. Lir > 10^^'^ Lq, corresponding to the bolometric luminosity of QSOs with Mb < —25) to be 
able to determine whether the luminosity function for ULIGs continues to evolve more strongly 
with increasing luminosity as does the luminosity function for QSOs. Deeper infrared surveys, 
such as are now being carried out with the Infrared Space Observatory (ISO) and in the future 
with the Wide-field Infrared Explorer (WIRE) and the Space Infrared Telescope Facility (SIRTF), 
when combined, should be able to measure the luminosity function of ULIGs out to the highest 
redshifts observed for QSOs, and thus will provide a clearer answer as to how well ULIGs mimic 
the luminosity function of QSOs. 



5.4. Clustering 

We have also investigated the possibility of clustering of ULIGs using the two-point correlation 
function ^(r) (Peebles 1980). At a comoving separation r, the two-point correlation function is 
given by 

= iVgg/ATg, - 1, (10) 
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where iVgg and ATgj. are galaxy-galaxy and galaxy-random pair respectively at a given comoving 
separation of r — Ar/2 < r < r + Ar/2. The errors were estimated from Poisson statistics 



In this calculation, the total number of a random sample which produces a minimum statistical 
noise was found to be ~ 100 x the total number of our sample. The results of the two-point 
correlation function analysis are plotted in Figure 6. As seen in Figure 6, most of the data points 
are consistent with no clustering except for one marginally consistent data point at r = 22 h"^ 
Mpc {h = Ho/100) where the correlation is 1.6 ± 1.2. 

Recently, Mo & Fang (1993) have detected strong clustering for quasars at 10 < r < 50 
Mpc. They have used about 750 quasars with mean redshift < z >= 1.5. To compare their result 
with ours, we have evolved Mo & Fang's result back from their mean redshift to our mean redshift 
< z >= 0.145 using the following three evolution models for the correlation function as proposed 
by lovino & Shaver (1988): 

i) Comoving model - The correlation amplitude A, defined by the expression ^(r) = Ar~^'^, 
has no dependence on redshift and the correlation function remains constant at all the time. In 
this case the clusters expand as the same rate as the Hubble expansion such as might be expected 
in biased theories of clustering (e.g. Davis et al. 1985). 

ii) Stable model - The correlation amplitude A has a (1 -|- z)^ dependence due to gravitational 
growth, so the correlation function has the form ^(r) oc (1 -|- z)~^-^. In this model the clusters keep 
the same proper distance with respect to the Hubble expansion. 

iii) Collapsing model - clusters collapse as fast as the Hubble expansion and the correlation 
amplitude A has a redshift dependence of < (1 -|- z)~^'^. 

For these calculations we have used A = 7A from the results of Mo &: Fang (1993). The 
resulting correlation functions for each evolution model are plotted in Figure 6. The collapsing 
model gives the poorest fit, and the comoving model (dashed lines) fits better than the stable 
model. This is consistent with Boyle &: Mo (1993) who find a similar clustering scale length for 
the low redshift quasars (z < 0.2). If the correlation function evolves according to the comoving 
model, then the above result gives some marginal evidence that both ULIGs and quasars may be 
derived from the same parent population. 



The following summarizes the main results from our IRAS 1-Jy survey of ULIGs selected 
from the IRAS FSC: 

1) A complete flux-limited sample of 118 ULIGs with /so > 1 Jy have been identified in a 
region of the sky at \b\ > 30°, S > —40° (~ 16,300 sq deg). The maximum observed infrared 
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luminosity is (Lir)max = lO^^-^'^ Lq, and the maximum observed redshift is z^ax = 0.268. 

2) The luminosity function, $ [Mpc~^ mag~^], for ULIGs can be approximated by a power law 
with spectral index a = —2.35 ± 0.30 over the decade luminosity range Lir = 10^^ — 10^^ Lq. In 

the "local Universe" (z ^ 0.1) the space density of ULIGs appears to be a factor of 1.5-2 times 
larger than the space density of optically selected QSOs at comparable bolomctric luminosities. 
However, if one accepts recent reported evidence for an increase in the local space density of 
QSOs, the space density of QSOs may be equal to, or even slightly exceed that of ULIGs in the 
local Universe. 

3) A maximum likelihood test suggests strong evolution in our sample. Assuming pure density 
evolution, $(z) oc (1 + we find a = 7.6 ± 3.2 . 

4) Examination of the two-point correlation function for our sample shows only marginal 
evidence for large-scale clustering, and only on size scales r ~ 22 Mpc. 
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FIGURE CAPTIONS 

Fig. 1. — AitofF projection in Galactic coordinates of all sources in the 1-Jy survey. The solid line 
represents the boundary of the survey area (see text). 

Fig. 2. — Finder charts for ULIGs in the 1-Jy Survey. The IRAS position uncertainty ellipse (iSa) 
in shown in overlay on the digital POSS red image (2'x2' field of view). A marker is centered on 
the identified ULIG. In those six cases with two marked sources, deeper optical and near infrared 
CCD images (Paper III) show that both galaxies appear to be strongly interacting (see text). 

Fig. 3. — Redshift and luminosity distributions of ULIGs in the 1-Jy survey. 

Fig. 4. — Mean spectral energy distributions of ULIGs (normalized at A = 60/Ltm). Arrows represent 
upper limits (2 a). 

Fig. 5. — Comparison of luminosity functions for different classes of extragalactic objects. 
References: IRAS 1-Jy Survey (this paper), IRAS RBGS (Sanders et al. 1998), Palomar-Green 
QSOs (Schmidt &; Green 1983), Markarian starbursts and Seyferts (Huchra 1977), and normal 
galaxies and cDs (Schechter 1976). The thick solid curve is a fit to the IRAS RBGS data points. 
The thick dashed line connects the IRAS 1-Jy data points. 

Fig. 6. — Two-point correlation functions of ULIGs. The comoving, stable and collapse models 
evolved back to z = 0.145 are plotted as dashed, long-dashed and dot-long-dashed lines respectively 
(see text). 



This figure "kim_pl_figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 
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TABLE 1— Continued 
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1 


.99 


la32 


0. 


242 


12 


.57 


12 


.63 


16333 + 4630 


16 


33 


22 


6 


46 


30 


55 


16; 


:34 


;50, 


9 


46:24 


;48 


72 


.30 


42 


.43 


0, 


.06 


0, 


.06 


1 


.19 


2 


.09 


la32 


0. 


191 


12 


.30 


12 


.35 


16468+5200 


16: 


:46: 


:49, 


,6 


52: 


:00: 


:52 


16; 


:48: 


;02. 


2 


51:55 


;40 


79 


.36 


39, 


.82 


0, 


.06 


0, 


.10 


1 


.01 


1, 


.04 


1232 


0. 


150 


11, 


.95 


12. 


.02 


16474+3430 


16: 


;47: 


;24, 


,6 


34: 


;30: 


;21 


16; 


:49: 


;14, 


,7 


34:25: 


:13 


56, 


.63 


39 


.12 


0, 


.13 


0, 


.20 


2, 


.27 


2, 


.88 


1332 


0. 


111 


12, 


.02 


12. 


.11 


16487+5447 


16: 


;48: 


;44, 


.1 


64: 


;47; 


;38 


16; 


:49: 


;47. 


,8 


54:42: 


:34 


82, 


.90 


39 


.23 


0, 


.07 


0. 


.20 


2, 


.88 


3 


.07 


1332 


0. 


104 


12, 


.07 


12. 


.12 


17028+5817 


17 


;02: 


;62, 


.7 


68: 


;17 


;63 


17; 


03; 


;41. 


,8 


58:13: 


;48 


87, 


.01 


36 


.86 


0, 


.06 


0, 


.10 


2, 


.43 


3, 


.91 


1332 


0. 


106 


12, 


.06 


12. 


.10 


17044+6720 


17 


;04: 


;29, 


.8 


67: 


;20: 


;34 


17; 


:04: 


;28. 


,5 


67:16: 


;34 


97, 


.99 


36, 


.12 


0, 


.07 


0, 


.36 


1 


.28 


0, 


.98 


1332 


0. 


135 


11, 


.97 


12. 


.13 


17068+4027 


17: 


:06: 


:53, 


,6 


40: 


:27: 


:12 


17; 


08 


;32, 





40:23 


;25 


64 


.76 


36, 


.14 


0, 


.08 





12 


1 


.33 


1 


.41 


1232 


0. 


179 


12 


.23 


12 


.30 


17179+5444 


17: 


:17: 


:54, 


,2 


54: 


:44: 


:51 


17; 


18 


;55, 


1 


54:41 


;50 


82 


.46 


35 


.04 


0, 


.08 





.20 


1 


.36 


1 


.91 


1332 


0. 


147 


12 


.07 


12 


.20 


20414-1651 


20: 


:41: 


:28, 


.4 


-16: 


:51: 


:08 


20; 


:44: 


;17. 


,4 


-16:40: 


;14 


29, 


.62 


-32, 


.27 


0, 


.66 


0, 


.35 


4, 


.36 


5, 


.25 


1232 


0. 


,086 


12, 


.06 


12. 


.14 


21208-0519 


21: 


;20: 


;61, 


.0 


-06: 


;19: 


;63 


21; 


:23: 


;28. 


,7 


-05:06: 


:59 


47, 


.25 


-36. 


.91 


0, 


.09 


0. 


.16 


1 


.17 


1, 


.66 


1132 


0. 


130 


11, 


.89 


12. 


.01 


21219-1757 


21: 


:21: 


:66, 


,1 


-17: 


:67 


:37 


21; 


24: 


;42. 


,5 


-17:44: 


;40 


32, 


.78 


-41, 


.66 


0, 


.21 


0, 


.46 


1 


.07 


1, 


.18 


3332 


0. 


112 


11, 


.69 


12. 


.06 


21329-2346 


21 


:32 


:54, 


2 


-23 


:46 


:03 


21: 


35 


;45, 


,4 


-23:32: 


;36 


26, 


.26 


-45 


91 


0, 


.05 


0, 


.12 


1 


.65 


2 


.22 


aa32 


0. 


125 


11 


.99 


12 


.09 


21477+0502 


21 


47: 


47, 


1 


05 


()2 


11 


21: 


50 


;17, 


6 


05:16: 


;14 


62 


.48 


-35, 


.57 





.09 





.16 


1 


.14 


1 


.46 


la32 


0. 


171 


12 


.12 


12 


.24 


22088-1831 


22: 


:08: 


:51, 


,5 


-18: 


:31: 


:50 


22; 


11 


;35. 


,5 


-18:17: 


;00 


37, 


.87 


-52, 


.26 


0, 


.09 


0, 


.07 


1 


.73 


1, 


.73 


la32 


0. 


170 


12, 


.29 


12. 


.31 


22206-2715 


22: 


;20: 


;41, 


.6 


-27: 


;16: 


;12 


22; 


23; 


;29. 


,5 


-27:00: 


;00 


24, 


.49 


-67, 


.19 


0, 


.10 


0, 


.16 


1, 


.75 


2, 


.33 


1132 


0. 


132 


12, 


.06 


12. 


.19 


22491-1808 


22: 


:49: 


:08, 


,7 


-18: 


:08: 


:24 


22; 


:51: 


;49. 


,0 


-17:52: 


;27 


45, 


.18 


-60, 


.96 


0, 


.06 


0, 


.66 


5, 


.44 


4, 


.46 


a232 


0. 


,076 


12, 


.06 


12. 


.09 


22541+0833 


22 


:54: 


:10 


,4 


08 


:33 


:12 


22; 


56 


;41, 


,3 


08:49: 


;15 


81 


.20 


-44, 


.55 


0, 


.09 


0, 


.18 


1 


.20 


1 


.48 


1132 


0. 


166 


12 


.11 


12 


.23 


23060+0505 


23 


()6 


:01 


8 


05 


()5 


:13 


23; 


08 


;34, 


2 


05:21 


;29 


81 


.65 


-49, 


.09 





.20 


0, 


.43 


1 


.15 





83 


3332 


0. 


173 


12 


.16 


12 


.44 


23129+2548 


23 


:12 


54, 


3 


25 


48 


()9 


23: 


15 


;21, 


6 


26:04 


:31 


97 


.40 


-31 


.99 





.08 


0, 


()8 


1 


.81 


1 


64 


la32 


0. 


179 


12 


36 


12 


38 


23233+2817 


23: 


:23: 


:20 


,5 


28 


:17: 


:49 


23: 


25 


;48, 


7 


28:34 


;19 


101 


.06 


-30, 


.64 


0, 


.13 


0, 


.28 


1 


.26 


2 


.11 


1332 


0. 


114 


11 


.83 


12 


.00 


23234+0946 


23: 


:23: 


:24, 


.0 


09: 


:46: 


:22 


23; 


:25: 


;56. 


,0 


10:02: 


;52 


90, 


.89 


-47, 


.43 


0, 


.06 


0, 


.08 


1 


.56 


2, 


.11 


la32 


0. 


128 


11, 


.99 


12. 


.05 


23327+2913 


23: 


;32: 


;43, 


,0 


29: 


;13: 


;29 


23; 


:35: 


;12. 


,5 


29:30: 


:05 


103, 


.67 


-30, 


.60 


0, 


.06 


0, 


.22 


2, 


.10 


2, 


.81 


1232 


0. 


,107 


11, 


.96 


12. 


.06 


23389+0300 


23: 


:38: 


:57, 


,7 


03: 


:00: 


:52 


23; 


:41 


;31, 


1 


03:17 


;31 


91 


.22 


-55, 


.22 


0, 


.09 


0, 


.35 


1 


.23 


1 


.17 


1132 


0. 


145 


11 


.99 


12 


.09 


23498 + 242:! 




49 




8 


24 




3') 


2:i 




;2(i 


4 


24:40 


:2() 


106 






28 





111 


n 


12 


1 


.02 


1 


45 


la 32 





212 


12 


29 


12 


41) 



Col 1: The IRAS object name in the Faint Source Database (FSDB). The prefix 'FSC is the standard designator for sources in the Faint Source Catalog. The prefix 
'Z' is used for objects that are not in the FSC — these objects are contained in the Faint Source Reject File (see Moshir et al. 1992). 

Cols 2-5: Right Ascension (hh:mm:ss.s) and Declination (dd: : ) of the IRAS source position as listed in the FSDB. Positions are given for epochs B1950.0 and 

J2aaa.a . 

Cols 6-7: Galactic latitude and longitude. The distribution of all 111 sources in Galactic coordinates is given in Figure 1. 
Cols 8-11: IRAS flux densities at 12/j,m, 25/j,m, GO/^m, and lOO^ra as listed in the IRAS FSDB. 

Col 12: Flux quality in each of the four IRAS bands. Numbers 1—3 arc the quality flags adopted for the IRAS catalogs {1 — upper limit, 2 — moderate quality, 3 — 
high quality). The letter 'a' means that the flux density was estimated by us after coadding all of the IRAS data available for the source using the ADDSCAN/SCANPI 
procedure (Helou et al. 1988); these 'low quality' measurements are typically 2— 3 detections. 

Col 13: Redshifts determined from our high resolution optical spectra (see §4.2). These redshifts are compared in the Appendix with our previous redshift measurements 
(e.g. Kim 1995) and with redshift data from the literature. 

Col 14: Far-infrared luminosity determined using the 60/J.m and 100/im IRAS fiuxes acording to the prescription given in Appendix B of Cataloged Galaxies and 
Quasars Detected in the IRAS Survey (1989) (see also Sanders & Mirabel 1996). 

Col 15: Infrared luminosity computed using the fiux in all four IRAS bands according to the expressions given in equations 1—3. 
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TABLE 2 
Redshifts 



Name 


our work 


previous values 




IRAS FSC 


•S^opt 




Ref 


Other Names 


00091-0738 


0.118 


0.118 


(18) 




00188-0856 


0.128 


0.129 


(18) 




00397—1312 


0.261 


0.262 


(19) 




00456—2904 


0.110 


0.110 


(1) 




r\r\AO(^ 0T01 

UU482— z7zl 


A 1 OA 

0.130 








(JIL) 1)4— 22*57 


Alio 

0.118 


0.118 


(1) 




n 1 1 AO A A 

(Jiibb— (J844 


Alio 

I). 118 


0.118 


(19) 




Uiiyy— z*>U7 


A 1 T^G 


0.156 


(1) 




Ulz98— U744 


A 1 0£? 

0.136 


0.136 


(18) 




01355—1814 


A 1 AO 

0.192 


0.191 


(19) 




01494-1845 


0.158 


0.157 


(1) 




01569-2939 


0.141 


0.140 


(1) 


TXS0156-296 


01572+0009 


A 1 (2 O 

0.163 


0.163 


(10) 


A 1 A1 /I T* A1 i AA1 

Mrkl(Jl4, FCjU 157+1)1)1 


02021—2103 


0.116 


0.116 


(1) 




02411+0353 


A 1 /i /i 

U.144 


0.144 


(19) 




02480—3745 


A 1 CT 

0.165 


0.165 


(19) 




03209—0806 


0.166 








03250+1606 


A 1 OA 

U.lz9 


0.129 


(19) 




ZjUoOzi+UUzo 


u.ioz 


0.152 


(18) 




A/inT/l OOA1 

04074—2801 


A 1 cr o 

0.153 








04103-2838 


0.118 


0.117 


(19) 




04313-1649 


0.268 








OtTT OA/11 

05020—2941 


0.154 








05024—1941 


0.192 








A CT 1 CT O AO /I 

05156—3024 


A 1 1^1 

0.171 








AC10A OCO/1 

(J5189— z5z4 


A A /I O 

0.042 


0.043 


(13) 




Ai^CAA 1 /^c:ao 

U7599+05U8 


A 1 /I A 

(J. 149 


0.149, 


0.148 (14), (19) 


Ptt'doaIatca 1 

[HB89j 0759+65 


AOOAI 1 OOA1 


A 1 GO 








AO /IT/I 1 1 O 1 O 

08474+1813 


A 1 /I C 

0.145 


0.146 


(19) 




08559+1053 


A 1 /I O 

0.148 


0.149 


(8) 




08572+3915 


0.058 


0.058 


(13) 




08591+5248 


0.158 








AAAOA 1 ACTAO 

09039+0503 


A IOC 

0.125 


0.125 


(19) 




09116+0334 


0.146 










U.iOU 


0.155 


(19) 




09539+0857 


0.129 








10035+2740 


0.165 








10091+4704 


0.246 


0.245 


(19) 




10190+1322 


0.077 


0.077 


(18) 




10378+1108 


0.136 


0.136 


(18) 




10485-1447 


0.134 


0.133 


(3) 




10494+4424 


0.092 


0.092 


(18) 




10594+3818 


0.158 


0.157 


(19) 
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TABLE 2— Continued 



Name 


our work 


previous values 




IRAS FSC 


Zopt 


-^opt 


Ref 


Other Names 


11028+3130 


0.199 








11095-0238 


0.106 


0.106 


(18) 




11119+3257 


0.189 


0.187 


(3) 




11130-2659 


0.136 






AMI 113-270 


11180+1623 


0.166 








11223-1244 


0.199 








11387+4116 


0.149 








11506+1331 


0.127 


0.127 


(18) 




11582+3020 


0.223 








Z11598-0112 


0.151 


0.151 


(18) 


[MHH96] J 1 20226-0 1 2948 


12018+1941 


0.168 


0.168, 0.169 


(4), (19) 




12032+1707 


0.217 


0.168 


(19) 




12072-0444 


0.129 


0.128 


(14) 




12112+0305 


0.073 


0.073, 0.072 


(13), (18) 




12127-1412 


0.133 








12265+0219 




0.159, 0.158 


(15), (18) 


3C273, PG1226+023, H1226+023 


12359-0725 


0.138 








12447+3721 


0.158 








12540+5708 


0.042 


0.042 


(3) 


UGC08058, Mrk231, VIIZw490 


13106-0922 


0.174 


0.174 


(19) 




13218+0552 


0.205 


0.190, 0.205 


(9), (19) 


[HB89] 1321+058, H1321+058 


13305-1739 


0.148 


0.148 


(1) 




13335-2612 


0.125 


0.125 


(1) 




13342+3932 


0.179 


0.179 


(17) 




13428+5608 


0.037 


0.037 


(13) 


UGC8696, Mrk273, IZw71, W851 


13443+0802 


0.135 


0.135 


(19) 




13451+1232 


0.122 


0.122 


(6) 


4C12.50 


13454-2956 


0.129 


0.129, 0.076 


(1), (18) 




13469+5833 


0.158 


0.158 


(19) 




13509+0442 


0.136 


0.189 


(19) 




13539+2920 


0.108 


0.108 


(18) 




14053-1958 


0.161 


0.162 


(1) 




14060+2919 


0.117 


0.117 


(17) 




14070+0525 


0.265 


0.265 


(19) 




14121-0126 


0.151 


0.150 


(19) 




14197+0813 


0.131 








14202+2615 


0.159 


0.158 


(17) 




14252-1550 


0.149 


0.150 


(19) 




14348-1447 


0.083 


0.082 


(13) 




14394+5332 


0.105 


0.105 


(17) 




14485-2434 


0.148 








15001+1433 


0.162 


0.163 


(18) 




15043+5754 


0.151 
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TABLE 2— Continued 





onr worlc 

\J LAJ. VV yjl. -LV 






IRAS FSC 


'^opt 


z t Ref 


Ot.Vipr IVpimp^ 


15130-1958 


0.109 


0.109 (18) 




15206+3342 


0.125 


0.125, 0.124 (14), (19) 




15225+2350 


0.139 


0.138 (19) 




15327+2340 


0.018 


0.018 (12) 


UGC09913, IC4553/4, VV540, Arp220 


15462—0450 


0.100 


0.101 (18) 




16090-0139 


0.134 


0.134 (18) 




16156+0146 


0.132 






16300+1558 


0.242 


0.242 (19) 




16333+4630 


0.191 


0.191 (19) 




16468+5200 


0.150 






16474+3430 


0.111 


0.111 (18) 




16487+5447 


0.104 


0.104 (18) 




17028+5817 


0.106 


0.106 (18) 




17044+6720 


0.135 


0.135 (19) 




17068+4027 


0.179 


0.072 (19) 




17179+5444 


0.147 


0.147 (19) 


TXS1717+547 


20414-1651 


0.086 


0.087 (18) 




21208-0519 


0.130 






21219-1757 


0.112 


0.113 (8) 


[HB89]2121-179 


21329-2346 


0.125 


0.125 (19) 




21477+0502 


0.171 






22088-1831 


0.170 


0.170 (19) 




22206-2715 


0.132 


0.131 (19) 




22491-1808 


0.076 


0.076, 0.078 (16), (18) 




22541+0833 


0.166 


0.166 (19) 




23060+0505 


0.184 


0.174 (7) 




23129+2548 


0.179 






23233+2817 


0.114 






23234+0946 


0.128 


0.128 (19) 




23327+2913 


0.107 


0.107 (18) 




23389+0300 


0.145 


0.145 (19) 


4C03.60 


23498+2423 


0.212 







References. — Redshifts: 1 - Allen et al. (1991), 2 - Arakelyan et al. (1971), 3 - Armus, Heckman 
& Miley (1989), 4 - Baan (1989), 5 - Bottinelli et al. (1987) 6 - Grandi (1977), 7 - Hill, Wynn- 
Willianis & Bccklin (1987), 8 - Hill ct al. (1988), 9 - Low ct al. (1988), 10 - MacKcnty & Stockton 
(1984), 11 - Melnick & Mirabel (1990), 12 - Mirabel (1982), 13 - Sanders et al. (1988a), 14 - 
Sanders et al. (1988b), 15 - Schmidt (1963), 16 - Soifer et al. (1987), 17 - Strauss & Huchra 
(1988), 18 - Strauss et al. (1992), 19 - Fisher et al. (1995) 
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TABLE 3 
Luminosity Function 



log Lir $ {V/V^)- 

[Lq] [Mpc-3 (Mag/2)-^] 

12.0-12.2 1.3 ±0.2x10"^ 0.55 ± 0.03 (63) 

12.2-12.4 3.9 ±0.7x10-8 0.55 ± 0.06 (35) 

12.4-12.6 9.6 ± 2.6 x 10"^ 0.50 ± 0.08 (14) 

12.6-12.8 2.2 ±1.0x10-9 0.61 ±0.15 (5) 

12.8-13.0 1.7 ± 1.7 X 10-1° 0.40 ±0.29(1) 



'^Number in parenthesis represents number of galaxies used for statistics. 
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This figure "kim_pl_fig2a.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig2b.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig2c.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig2d.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig2e.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



This figure "kim_pl_fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/astro-ph/9806148vl 



